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ABSTRACT 

The Lagrange planetary equations are used to study to secular evolution 
of a small, eccentric satellite that orbits within a narrow gap in a broad, self- 
gravitating planetary ring. These equations show that the satellite's secular 
perturbations of the ring will excite a very long-wavelength spiral density wave 
that propagates away from the gap's outer edge. The amplitude of these waves, 
as well as their dispersion relation, are derived here. That dispersion relation 
reveals that a planetary ring can sustain two types of density waves: long waves 
that, in Saturn's A ring, would have wavelengths of A ~ (9(100) km, and short 
waves that tend to be very nonlinear and are expected to quickly damp. The 
excitation of these waves also transports angular momentum from the ring to the 
satellite in a way that damps the satellite's eccentricity e, which also tends to 
reduce the amplitude of subsequent waves. The rate of eccentricity damping due 
to this wave action is then compared to the rates at which the satellite's Lindblad 
and corotation resonances alter the satellite's e. These results are then applied 
to the gap-embedded Saturnian satellites Pan and Daphnis, and the long-term 
stability of their eccentricities is assessed. 

Subject headings: planets: rings 
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Introduction 



The following considers the secular gravitational perturbations that are exerted between 
a small eccentric satell ite and a nea rby planetary ring. This investigation is also a followup to 
the study described in lHahnl (120071 ). which examined the secular evolution of a small inclined 
satellite. There it was shown that the inclined satellite would launch a very long-wavelength 
spiral bending wave at the ring's nearby edge. Since the excitation of the bending wave also 
communicates in-plane angular momentum from the satellite to the ring, that ring-satellite 
interaction resulted in a very vi gorous damping o f the satellite's inclination. A related 
problem was also considered by iGoldreich fc Saril (120031 ). who showed that an eccentric 
perturber orbiting in a gap in a pressure-supported gas disk will launch density waves at 
the gap's edge having such a long wavelength that a global standing wave emerges in the 
disk. The work described below considers a related but distinct problem, that of an eccentric 
satellite orbiting in a narrow gap in a self-gravitating planetary ring. Here we show that 
the secular perturbations from the eccentric satellite also launches relatively long-wavelength 
density waves at the gap's outer edge. The principal goal of this study will be to derive the 
amplitude of these waves and their dispersion relation, which in turn will yield other useful 
properties, such as the wavenumber and the waves' group velocity. Since the mathematics 
of this density wave problem is very similar to that already developed for bending waves, 
the derivations presented here will be succinct. How ever, the rea der interested in a more 
verbose description of a similar problem is referred to iHahnl ( 120071 ). 



Another goal of this effort will be to derive the rate at which the satellite's eccentricity 
e is altered by the excitation of these density waves at the gap's edge. This is of interest 
because the satellite's e also varies due to its many Lindblad resonances in the ring, which 
tend to pump up the s atellite's e, while its many corotat ion resonances in the ring tend to 
damp the satellite's e ( IGoldreich fc Tremaindll98ll . Il982f ) . Although e-damping due to the 
corotation torque can dominate over e-excitation from the Lindblad torque, the long-term 
stability of the satellite's eccentricity is still uncertain since the corotation torque is operative 
only when these particular resonances are not saturated. Our purpose here is to determine 
whether the secular interaction described below is also a stabilizing process that leads to a 
net damping of the satellite's eccentricity. 

The following Section begins with the Lagrange planetary equations, which are used to 
derive the amplitude of these waves and their dispersion relation. Section [3] then uses those 
results to determine the rate at which the satellite's eccentricity evolves in response to the 
wave it launches at the gap edge. Those analytic results are then confirmed in Section @] via 
a numerical simulation of these waves, with conclusions reported in Section [51 
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2. Equations of motion 

Consider a planetary ring that is perturbed by a single satellite, with both orbiting an 
oblate planet. The Lagrange planetary equations give the rates at which a ring particle's 
orbital eccentricity e and longitude of periapse u vary with time t due to these perturbations 



(IBrouwer fc Clemencelll96ll ; iMurray fc Dermottlll999f ): 



dR ■ 1 dR 

and lu ~ — - — — — , (1) 



na 2 e dui na?e de 

where R is the disturbing function for a small ring particle having a semimajor axis a and 
mean motion n ~ y^GM/a 3 , where G is the gravitation constant and M is the mass of the 
central planet, and all eccentricities are small, e< 1. The total disturbing function for a 
ring particle is R = -Rdisk + -Rsat + -Robh where the three terms account for the gravitational 
perturbations that are due to the ring's gravity (which is treated here as a broad disk), the 
satellite's perturbations, and that due to the planet's oblate figure. The particle's equation 
of motion is thus the sum of three parts: 

e= e| disk +e| sat and Q = i>| disk + + £>| oM , (2) 

noting that oblateness does not alter eccentricities. And because we are only dealing 
with the system's secular p erturbations, the system's semimajor axes a are all constant 



( IBrouwer fc Clemencdll96ll ). 



The amplitude of a spiral density wave that is in a steady-state does not vary with time, 
so the disk eccentricities obey e(a) = throughout the disk. A persistent spiral pattern must 
also rotate with a constant angular velocity u, so 

^ I disk = _ ^Isat ( 3a ) 

and uj = Co ,. . + uj\ . + Cb\ = constant. (3b) 

I disk I sat I obi v ' 

These equations are used to obtain the wave amplitude e(a) throughout the disk, and the 
waves' dispersion relation u(k) as a function of the wavenumber k. 



2.1. wave amplitude 

First calculate the rate at which the planetary ring perturbs itself. This ring is treated 
as a broad disk composed of many narrow, concentric annuli that have a mass 5m(a), ec- 
centricity e(a), and longitude of periapse uj(a) that are regarded as functions of the rings' 
semimajor axes a. Each annulus also has a half-thickness h that is due to the ring particles' 
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dispersion velocity. Now suppose that the annulus at a is perturbed by another annulus of 
mass 5m' and radius a'; the perturbed annulus will then have a disturbing function 



5R 



G5m' 
4a 



-f(a)e 2 + g(a)ee cos(u — Co' 



(4) 



where a, e, co> are the orbit elements of the perturbed annulus, the primed quantities refer to 
the perturbing annulus, and the f{a) and g(a) functions are 



/(«) 
and g(a) 



(a) - 6f) V6$,(a) 
-a6^(a) + 6f)V^ 1 / ) 2 (a), 



(5a) 
(5b) 



which depend on the semimajor axis ratio a = a' /a and the disk's dimensionless scale 
height f) = h/a <C 1, which is presumed small. This disturbing function is derived in lHahn 
(120031 ). and it differs somewhat from the more familiar disturbing function for a point-mass 
perturber due to the annuli's finite thickness h, which also softens the Laplace coefficients 
that appear in 5R. Those softened Laplace coefficients are 



a) 



cos(m<p)d(p 



[(l + a 2 )(l + f) 2 )-2acos<^' 



(6) 



Note that when the annuli are infinitesimally thin, t) — > 0, and the softened Laplace co- 



efficients b 



s 



oc, 



are equivalent to the conventional Laplace coefficients b. 



(m). 



a: 



for which 



f(a) — > a&g^a) an d g>(a) — >■ —abf) 2 (a), and the di sturbing function SR becom es equivalent 



to that due to a point-perturber of mass 8m' (e.g., iBrouwer &: Clemencdll96ll ). 



The perturbing ring's mass is 5m' = 27T<y'a'da' where a' = cr(a') is the mass surface 
density of the perturbing annulus of radius a' and radial width da'. The perturbing ring's 
semimajor axis is written a' = a(l+x'), where x' = (a'— a) /a = a — I is the fractional distance 
between the perturbing ring a' and the perturbed ring a. Also define the normalized disk 
mass as Hd{°) = vrcra 2 /M, so that the factor G5m' /4a in Eqn. P| becomes n' d (na) 2 dx' /2a 
where aV = da' /a is the perturbing ring's fractional width. The disturbing function for ring 
a due to perturbations from ring a', Eqn. (j3J, then becomes 



-f(x')e 2 + g(x')ee' cos(u — Co') 



dx' 



(7) 



where f(x') and g(x') is shorthand for Eqns. (jSJ) evaluated at a = 1 + x'. Inserting this into 
Eqn. then shows that ring a' alters the eccentricity of ring a at the rate 



5e 



1 d(SR) 
na 2 e duo 



-^' d na 1 g(x')e sin(u; — uj')dx . 
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2.1.1. ring-disk evolution 

The total rate at which the entire disk alters the eccentricity of ring a is the above with 
x' integrated across the disk, so e| disk = f disk 5e. For the moment, consider a one-sided disk, 
one that orbits wholly exterior to the satellite, where A is the fractional distance between 
the satellite's orbit and the disk's inner edge. The geometry is sketched in Fig. [U which 
shows that the integration variable x' then ranges from — x to +00, so the ring's eccentricity 
varies at the rate 

e|disk=/ <te = -n/ rf d (x')a- x g{xys\T\{u)-u')dx'. (9) 

./disk ^ J-x 

This integral will be dominated by the contributions from nearby annuli that lie a small 
distance x' away. In the \x'\ <C 1, f) <S 1 limit, the softened Laplace coefficients that appear 
in the g(a) function, Eqn. (j5bj) . are 



6 »M = ^TW) (10a) 

4 



and ^ * 3rr(x' 2 + Wf (1 ° b) 



(lHahnll2003l ). so a 1 g(a) in Eqn. (j5j) becomes 



a 1 g(a) ~ ^— (11) 

in this approximation. Due to the steep dependence of g(a) on x', we can replace the eccen- 
tricity e'(x') and disk mass ^{x') in Eqn. ([9]) with their values evaluated at the perturbed 
ring, which lies at x' = 0, so e' ~ e and ~ fx d = 7iaa 2 /M, which are then pulled out of 
the integral so that 

1 r x l2 -2\) 2 . , . 

e ldisk - "-^en y ^ p - ^ 2)2 sin( W -u)dx. (12) 

A spiral wave will have a wavelength A ~ 2ir/\k\, where k = —du/da is the wavenumber 
of the spiral density wave. Thus the to — Co' in Eqn. (1121) can also be written as 

Co — Co (a) — — I k{r)dr. (13) 



However, most of the contribution to the integral in Eqn. (j!2p will be due to nearby annuli 
that lie a wavelength A away. But if k also varies slowly with a, then it can be treated as a 
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constant over that wavelength, so Eqn. (fT3l) is u — uj' ~ —k(a — a' 
then becomes 



kax', and Eqn. (fl~2l) 



I disk 



kaiid&n 



71 



H 2 



1 



\k\ax 



(y 2 + 2H 2 ) 2 



sin(y)dy = A' H (\k\ax)ka/j,den 



7T 



(14) 



after replacing the x' integration variable with y = \k\ax', which is the distance from the ring- 
edge in units of 2tt wavelengths, and replacing the disk's scale height \) with H = y/2t)\k\a, 
which is roughly the disk's vertical thickness in wavelength units. In the above, the function 
A\z) is 

(y 2 + H 2 ) 2 Sinydy ' (15) 

noting that Eqn. (1141) is also odd in y. We will be interested in a disk whose vertical 
thickness is small compared to the wavelength, so H <C 1, and A ' H (z) ~ sin(z)/z — Ci(z), 
where Ci(z) is the cosine integral of lAbramowitz Stegunl (119721) . Far dow nstream, where 
z » 1, G\{z) ~ sax{z)/z — cos(z)/z 2 + 0(z~ 3 ) (lAbramowitz fc Stegunl Il972l ). so 



AM 



cos(^) 



(16) 



downstream. Evidently, A' H (z) is very similar to the Ah(z) function of iHahnl (120071 ). which 
is also plotted in Fig. 2 there. 

Also keep in mind that Eqn. ([Uj) was derived with the understanding that the wavenum- 
ber k varies little over a single wavelength. Section [2.2.31 will quantify when that assumption 
breaks down. 



2.1.2. ring-satellite evolution, and the wave amplitude 



The ring at semimajor axis a is also perturbed by the satellite, and the disturbing 
function R s due to the satellite is Eqn. (j4j) with 5m' replaced by the satellite's mass m s , 
which can be written 



R.H 



1 



-u s [na) 



-f(a)e 2 + g(a)ee s cos(a; — u s 



;i7) 



where fi s = m s /M is the satellite's mass in units of the central planet's mass, and a = a s /a = 
(1 + A + x)^ 1 ~ 1 — (A + x). The satellite's perturbation causes the ring's eccentricity to 
vary at the rate 

1 BR, 1 , . 

18 



I sat 



na 2 e dQ 



-/j, s ng{a)e s sm[uj - u s ) 
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The perturbed ring lies a fractional distance x + A away from the satellite (see Fig. [I]), with 
both well separated so that A>^. Since the ring also lies in the wave-excitation zone near 
the satellite, 

~ 2 
g(a) ~ -abfj 2 (a) ~ --^-—^ (19) 



> + A) 



by Eqns. ( I5bl) and (llOaj) . Also write the longitude difference in Eqn. (TT5T) as — cD s ~ 
—kax + (p , where the angle <p allows for the possibility that the annulus nearest the satellite 
at x = has a longitude of periapse that differs from the satellite's longitude u) s by angle 
4> - Thus 

e I sat - 27r (^+ A)^ sin(A:aX ~ ^ (20) 
is the rate at which the satellite alters the ring's eccentricity. 

When the wave is in steady-state, the two e-excitation rates, Eqns. (I14p and (I20p . are 
balanced, which provides the amplitude of the density wave: 



e(z) _ \k\an s sin(z - ..,.,,„, ? 
e s ~ 2 ^ d {z+\k\a^fA' H (z) [ j 

where Sk = sgn(/c) and z — \k\ax is the downstream distance in units of 2tt wavelengths. Now 
recall that A' H (z) ~ cos(z)/z 2 far downstream where z ^> 1 (Eqn. [IB]) . Since the downstream 
wave amplitude e(z) should be a finite constant for all z ^> 1, this then tells us that the 
longitude offset is O = ±vr/2, so e/e s ~ — (|&|a/2)(/i s //id)sfc sin0 o . These eccentricities must 
also be positive, so sin0 o = ±1 = — Sk, and the density wave amplitude then becomes 

e \ k o\a^s ^ 



e s 2fi d 

where \ko\ is the initial wavenumber at x = 0, where the wave is excited at the disk's inner 
edge. An identical expression was also obtained in iHahnl (120071 ) for the amplitude of the 
bending wave that an inclined satellite would launch at the ring's edge. 



2.2. dispersion relation 



Further use of the wave amplitude, Eqn. (122!) . requires knowing the initial wavenumber 
ko, which is obtained from the waves' dispersion relation, Eqn. ( I3bl) . The first term in that 



equation, uo 



disk : 



is the rate at which the disk drives its own precession. The rate that a 



single annulus at a precesses due to the secular perturbations from another annulus at a' is 
1 d{8R) 1 



5uj 



na 2 e 



de 



-u d na 



f(x') +g(x' 



e'(x') 



cosfu; — U)') 



dx' 



(23) 
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where 5R is Eqn. (J7j). The total precession rate due to the disk's self-gravity is £)| disk = 
J k Su, where the integration proceeds across the entire disk. Again, the integrand is a 
steep function of x' due to the softened Laplace coefficients that are present in the / and 
g functions, which allows us to replace the quantities e'(x') and fi' d (x') with the constants e 
and Hd and to pull them out of the integral. And since /(a/) — —g{x') when \x'\ <C 1 (from 
Eqns. |5]and HO]), 



LO 



disk 



1 

-fjL d n 

7T 

2 

-fjL d n 

7T 



X 



n 



,,/2 



[.v- -2f) , 

,2 _ 2[) 2 



2fj 2 

— [1 — cos (| A; | ax')] da/ 



x ■ 



x' 2 - 2f) 2 ) 2 
B' H ( | k | ax) | k | a/irf/i 



sin 2 ( | k | ax' /2) dx' 



where 



B' 



sin 2 (y/ '2) dy. 



(y 2 + H 2 ) 2 

When the disk is much thinner than the wavelength, H = \^2l)\k\a 1, and 



B'h(z) 



1 
2 



-Si(^) 



cos z 



71 



TCZ 



(24a) 

(24b) 
(24c) 

(25) 
(26) 



where Si(z) is the sine integral of lAbramowitz &: Stegunl (19721 ). In this limit the B' H (z) 
function is identical to the Bh(z) function of iHahnl (120071 ). with both taking values of 
1/2 < B' H < 1. And far downstream, where z —>• oo, the B' H integral evaluates to 



B'h = lim B' H {z) 



(27) 



for any value of H, which, by the way, does differ a bit from the Bh function of lHahnl (120071 ). 



The function B' H (z) is a dimensionless measure of the rate at which the disk drives its own 
precession, and this quantity becomes small when H > 1, or when the disk thickness exceeds 
the wavelength. Evidently, a thick disk having H 3> 1 is less likely to sustain a density wave. 

The satellite's gravity also precesses the planetary ring, at the rate 

1 dH, 1 



LO 



I sat 



na 2 e de 

Us 



2ir(x + A) : 



+ 



f(a) + g(a)—cos(u - io s ) 



[id sin ( A; 


ax) 


7T 


k 


a(x + 


A) 2 



e 
n 



(28) 
(29) 



when e s /e is replaced with Eqn. (1221) . The disturbing function associated with the central 
planet's oblateness is also 



R, 



obi 



: J 2 e 2 



an 



(30) 
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wher e J2 is the planet's second zonal harmonic and R p is the planet's radius ([Murray fc Dermott 



19991 ). so precession due to oblateness is 



1 OR, 



to 



obi 



obi 



na 2 e de 



-,/:, I 



n 



l--(x + A) 



obi 



(31) 



where & s \ ohl = (3J2/2)(-R p /a s ) 2 n s is the rate at which the satellite's orbit precesses due to 
oblateness, with n s being the satellite's mean motion. Summing Eqns. (I24cl) . (129j) . and (13T1) 
then provides the dispersion relation for these spiral density waves, 



a>(|fc|) ~ D' (z) fi d \k\an + 



l- 7 -(x + A) 



I obi 



(32) 



which is the angular rate at which the spiral density pattern rotates. The D' function in the 
above is 



D'{z) = B' H (z) + 



sin z 



(33) 



ir(z+ \k\aA) 2 ' 

where D'(z) is again numerically similar to the D(z) function of iHahnl (120071 ). which is also 
plotted in Fig. 2 there, which shows that D'(z) takes numerical values of 1/2 < D(z) < 1 
when H <C 1. 



2.2.1. long and short density waves 

When these density waves have propagated far downstream of the satellite, the middle 
term in Eqn. (1321) . which is the rate at which the satellite precesses the disk, becomes even 
more negligible with distance from the satellite, and the D'(z) function becomes D' ~ B'^ = 
e~ H where H = y/2t)\k\a should now be regarded as a dimensionless wavenumber. Thus the 
downstream dispersion relation is 

u(\k\) ~ fi d \k\ane-^ kla + Co\ qU (x). (34) 

This dispersion relation can then be expressed in a more convenient dimensionless form via 
the combination 

"* = ^M\k\)-t>U = He- H - (35) 

Figure [2]plots this dimensionless dispersion relation uj* versus the dimensionless wavenum- 
ber H, which shows that uj* has a maximum value of c^ max = exp(— 1) ~ 0.368, which occurs 
for a wavenumber H = 1. This figure also shows that as long as the density wave's spiral 
pattern rotates at an angular rate uj* < c<J max , then the disk can sustain two types of waves: 
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long waves that have a wavenumber H < 1, and short waves that have H > 1, noting that 
these waves are called such since their wavelength is A = 2ir/\k\ = 2\/2irh/H. This upper 
limit on u* also tells us that the disk can sustain these density waves when the spiral density 
pattern does not rotate too fast, namely, when 



u < 



+ Cj 



obi 



(36) 



see also lHahnll2003l ) 



2.2.2. group velocity 



The waves' group velocity is c g = du/dk (lToomrelll969l ; IShulll984j ; iBinney fc Tremaine 

19871 ). which becomes 

c g = S kg^ = Sjfe(l - H)e~ H n d an (37) 

when the downstream dispersion relation is differentiated (Eqn. [51)1 ; this is the rate of the 
spiral wave's radial propagation. Since the satellite is launching outward-propagating density 
waves from the disk's inner edge, c g > 0, so Eqn. ([37]) implies that the satellite can launch 
long H < 1 waves that have Sk = sgn(k) = 1, or short H > 1 waves that have Sk = — 1. 
Spiral density waves having k > are called trailing waves, since the more distant parts of 
a spiral arm trail in longitude, while leading waves have k < 0. 

Section [223] will show that any waves excited by the Saturnian satellites Pan and Daph- 
nis, both of which inhabit narrow gaps in Saturn's main A ring, would have a dimensionless 
wavenumber H -C 1. This means that these satellites could launch l ong trailing waves that 



have Sk — 1 that would propagate outwards at a rate c 9 ~ /i^cm (e.g., lHahnll2003l ). And since 
sin0 o = —Sk = —1 (see Section [2. 1.2j) . then <p = —n/2, which means that the longitude of 
periapse at the disk's inner edge would trail the satellite's longitude by 90°. 

Later, Eqn. (|47p will also show that the wavenumber \k\ and hence H will increase as 
the wave propagates downstream. So it is possible that a wave might travel far enough 
for the wavenumber H to increase beyond unity. If that happens, then Eqn. (137)) suggests 
three possible outcomes: the trailing wave can continue to advance further downstream 
with c g > as a short (H > 1) leading (sk = —1) wave, or it can reflect (c g < 0) as a 
short (H > 1) trailing (s& = +1) wave, or it might do both by spawning both types of 
wavetrains. Note, though, that the short wave would propagate at a substantially slower 
rate, \c g \ ~ He~ H \i&an, since H > 1. This site in the disk where H = 1 and c g changes 
sign is a turning point for long waves. This site is also known as the Q-barrier, since its 
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location depends upon the value of Toomre's stability parameter Q ~ v^nj-nGo = t)/nd 
where Vd = hn is the particle disk's dispersion velocity (jToomrelll964j ). 



2.2.3. wavenumber and wavelength 

The wavenumber k is obtained after calculating the satellite's precession rate tij s . When 
the system is in steady state, both the satellite and the spiral wave precess at the same rate, 
Cj s = uj(\k\), which then provides a simple relation for k that depends only on the system's 
physical constants, namely, fid, A, and J2. 

The satellite's longitude of periapse Cb s precesses due to perturbations from the disk and 
the central planet oblateness, with this precession occurring at the rate 



10 



s I disk + ^lobl 



where Co. 



I disk = idisk^s * s ^ ne satellite's precession rate due to the entire disk, where 

1 



Suj, 



/ 9{<x 



a 



1 COS{U! s 

e.q 



Co') 



dx' 



(3J 



(39) 



is the satellite's precession rate due to a disk annulus of radius a' and mass 8m' . This is 
obtained from Eqn. ( |23j) when n, a,e, Co is replaced with n s ,a s ,e s ,Co s , and the separation 



x 



x' + A. The satellite's precession rate due to the entire disk is 



disk 



8Cb s ~ -fi d n 

disk 71 

ttA 2ttA 2 



X 



A) 



-2 



— cos( kax' 



dx' 



S(\k \aA) 



where 



S(\k \aA) = \k aA\ 



sin(y)dy 



(40a) 
(40b) 

(41) 



/ (y+\k \aA) 2 

is another dimensionless function that depends on the t he init i al wav enumber \ko\ and the 
gap width A. This quantity is also plotted in Fig. 2 of iHahnl (120071 ). which shows that S 
takes numerical values of < S(\ko\aA) < 1. 

The first term in Eqn. fl40b|) is the rate at which the undisturbed disk precesses the 
satellite's orbit, while the second term is the additional precession that the satellite experi- 
ences due to the density wave in the disk. Note, though, that if the satellite instead orbited 
at the center of a narrow gap in the disk, then the first term in Eqn. f!40bj) would be doubled 
due to the disk matter orbiting interior to the satellite. One would also expect additional 
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precession to occur due to any density waves that might be launched in this interior disk 
However, it will be shown below that this contribution is unimportant. With this in mind 
Eqn. f!40bj) is generalized to account for a possible inner disk by writing 

e^ d n s i fi s n s 



oj. 



disk 



+ 



;S(\k \aA) 



(42) 



ttA 2ttA 2 

where e = 1 if the disk lies entirely interior or exterior to the satellite, and e = 2 if the satellite 
instead orbits in the center of a gap whose fractional half-width is A. The satellite's total 
precession rate then becomes 

e\i d n s ( fi s n s 



10, 



disk 



+ 0J t 



obi 



ttA + 2ttA 2 



^d^olaA) + u s 



I obi ' 



(43) 



When the disk and satellite are in steady-state, the satellite and its spiral wave pattern 
both precess in concert, so uj s = u(\k\), which becomes 



/'.- 



where 



nD'(z)\k\aA = e + ^ (l + + ^f(\ k aA\, z) 
fji d \ A/ 

f(\k aA\,z)=S(\k \aA) 



2/i d A- 
\k aA\ 2 



(\k aA\+z) 2 

is another function that takes values of — 1 < / < 1, and 

2 



- / R p A 



(44) 



(45) 



(46) 



is called the critical disk mass (jHahnl 120071 ). 



Two additional assumptions will then provide a simple expression for the wavenumber 
k. First, assume that the right term in Eqn. (j4"4"|) is small compared to the middle term, 
which requires the satellite's mass to be sufficiently small, namely, that fi s 2fi c A. Second, 
assume that the dimensionless wavenumber obeys H <C 1, which means that D'(z) — > 1 
downstream. Note, though, that D'(z) is not unity in the wave laun ch zone; rat her it takes 
values of 1/2 < D'(z) < 1 over the initial wavelength (see Fig. 2 of iHahnl 120071 ). However, 
a sufficiently accurate result is obtained wh en D'(z) is replaced with its average value over 
the first wavelength, D ~ 0.87 (jHahnl 120071 ). These assumptions will be confirmed below in 
Section f2.2.4[ With these assumptions in hand, Eqn. fT44|) is then yields the wavenumber 



\k(x)\ 



nDaA 



( \ X 

A«d V A 



with 



l^o I 



nDaA 



(47) 



(48) 
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being the initial wavenumber evaluated at x — 0. 

The first wavelength A is obtained from J Q A ° \k\da = 2ir. When fi c -C fid, that integral 
yields 

A ~ 2n 2 DaA/e. (49) 

However, when fi c > fid, (which occurs when the central planet is too oblate or the disk mass 
is too small), the wavenumber k varies with distance x, which violates a key assumption in 
the derivation of the wave amplitude (Section |2.1.ip . The numerical experiment described 
in Section 14.31 will also show that the wave amplitude (Eqn. derived below) is smaller 
than expected in this case. Nonetheless, that experiment also shows that the formula for the 
expected wavenumber, Eqn. (I4TI) . is quite reliable when fi c > 



2.2.4- checking the assumptions 



2008). Pan orbits 



Porco et all 120051 ) 



These results will be applied to the Saturnian satellites Pan and Daphnis that orbit 
i n narrow gaps in Saturn's main A ring. Pan has a mass ix x = 8.7 x 10 ~ 12 Saturn masses 
(jPorco et al.ll2007l ) and a semimajor axis a s = 133, 584 km (IJacobson et a 
near the center of the Encke gap whose half-width is Aa = 162.5 km ( 
so its fractional half-width is A = Aa/a s = 0.0012. Saturn's radius is R p 
and its second zonal harmonic is J<i = 0.0163, so the critical disk mass is fi c 
Consequently, the quantity 2/z c A 
that fi s <C 2fi c A. 



60, 330 km 
7.9 x 10~ 8 . 

1.9 x 10~ 10 , which does indeed satisfy the requirement 



This same as sumption is also satisfied by Daphnis, whose mass is fx s = 1.5 x 10~ 13 
( IPorco et al.l 120071 ). but by a smaller margin. Daphnis orbits in the Keeler gap, which itself 
is about 2900 km be yond the Encke gap. The Keeler gap is also about 8 times narrower 
than the Encke gap (ITiscareno et al.l 120051 ). so A ~ 1.5 x 10~ 4 , which reduces the critical 
disk mass to fi c ~ 1.2 x 10~ 9 . Consequently, fi s < 2/i c A is still satisfied, but only by a factor 
of ~ 2.5. 

The other requirement is that H = \/2l)\k\a << 1 . The ring's vertical thickness is 
about h ~ 30 m near the Encke gap (ITiscareno et al.ll2007l ). so the ring's fractional thickness 
is f) = h/a s ~ 2 x 10~ 7 ; a comparabl e thic kness can also be inferred from the A-ring 
viscosity that is reported in IPorco et al.l (120071 ). The ring's surface density here is a ~ 46 
gm/cm 2 (ITiscareno et al.l 120071 ). so the normalized disk mass is fid = naa^/M ~ 4.5 x 10 -8 , 
which means that Pan's initial wavenumber is \ko\a ~ 1100 by Eqn. (1481) . while Daphnis 



has |fc |a ~ 5000. Consequently, H = y/2t)\k\a ~ 0(W 



3 to -4> 



<< 1 is very well satisfied. 



And since H < 1, this means that these satellites will launch long trailing density waves 
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that propagate radially outwards from the gap edge, provided these satellites have nonzero 
eccentricities. 



2.2.5. waves in an interior disk 



A satellite that orbits in the center of a gap in a broad planetary ring might also excite 
a disturbance in the disk material that orbits interior to the satellite. The dispersion relation 
for any density waves launched at the inner gap edge is 



\k{x)\ 



vrDalAI 



X 




A 


)] 



(50) 



which may be derived via the method that is described in footnote 1 of iHahnl (120071 ). This 
dispersion relation is identical to Eqn. (HTl) except for the sign on the term that is proportional 
to the critical mass. Since the right hand side of Eqn. (1501) must be positive wherever the 
wave propagates, this equation tells us that these density waves are only allowed inwards a 
distance \x\ < x- in , where 



IAI 



(51) 



is the distance of the waves' maximum inwards excursion. For Pan, this distance evaluates 
to x m ~ 0.1A, which corresponds to a physical distance of x- in a s ~ 15 km, which is only 
a tiny fraction of the wave's initial wavelength that is calculated below in Section 12.31 In 
short, the A ring material that orbits interior to Pan is unable to sustain this type of density 
wave. 

Note, however, that x- in is a bit larger for the ring material that orbi ts interior to Daphnis 
and the Keeler gap. H owever, Daphnis' small mass (IPorco et al.l 120071 ) and low eccentricity 
(IJacobson et al.l 120081 ) results in waves of such low amplitude that this issue is moot. 



2.3. wave amplitude, continued 

Plugging the initial wavenumber, Eqn. (T4"gj) . into Eqn. (I22p then yields an expression for 
the amplitude of the density wave that is now written only in terms of the system's known 
physical parameters: 

e s ~ 2nDfi d A ' 1 ' 

However, keep in mind that this derivation requires that the wavelength vary little over the 
first wavelength, which in turn requires \i c <C [id- The simulations described in Section 14.31 
will show that Eqn. (1521) overestimates the wave amplitude when \i c > fi d . 
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Pan's eccentricity is e s = 1.4 x 10~ 5 (IJacobson et al.ll2008l ). so Eqn. (152]) predicts eccen- 
tricities of e ~ 1.5 x 10~ 6 due to the density wave that Pan excites at the outer edge of the 
Encke gap, since e = 2, /i c = 7.9 x 10~ 8 , and fid ~ 4.5 x 1CT 8 (see Section [2.2. 4p . However, 
this is actually an overestimate, since fi c j fid — 1-8 is n °t small, which is a key assumption in 
the derivation of Eqn. (|52|) . In fact, the simulation described in Section T4.3I will show that 
Eqn. (152]) overestimates the amplitude of Pan's waves by a factor of 7 ~ 4, so the epicylic 
amplitude due to this wave at the outer Encke gap is only Ar = ea s /7 ~ 50m, which is likely 
too small to be seen by a spacecraft such as Cassini. And since this wave's initial wavenum- 
ber is \k \a ~ 1100 (Section |2.2.4[) . its initial wavelength would be Ao = 2n/\k \ ~ 760 km, 
which is also far longer than the wavelengths of any of the more familiar density waves that 
satellites routinely launch at their many mean-motion resonances in the rings. 

Recall that the wavenumber \k\ increases with distance x (see Eqn. I4T|) . so the wave- 
length A shrinks as the wave propagates outwards. Note that the A ring's outer edge is 
Aa ~ 3200 km away from the Encke gap, which corresponds to a fractional distance of x = 
Aa/a s ~ 0.024. Inserting this into Eqn. (1471) then yields a wavenumber of \k\a ~ 1.2 x 10 4 , 
which corresponds to a wavelength of Ao = 2ir/\ko\ — 70 km when the wave hits the outer 
edge of the A ring. This presumes that the wave was not damped en route by the ring's 
viscosity, but Section 12.3.11 will show that the viscous damping of these long outbound waves 
is modest. The reflection of these undamped waves near the outer edge of the A ring is also 
discussed in Section 



The fractiona l variation in the disk's surface density due to the density wave is (IBorderies et al 



19851 : iHahnl 120031 ) 



Aa d(ea) .„ dui . 

~ — - — cos(# — uj) + ea— sm(v — cu), (53) 

o"o Oa da 

where 6 is the longitude in the disk and <To is the surface density of the undisturbed disk. 

The first term is negligible since the downstream eccentricity e is constant, so the magnitude 

of the fractional variation in surface density variations due to the wave is dominated by the 

second term, which is 

ea\k\ (54) 



00 

since k = —du/da (see Eqn. [TBI . Since the wave launched by Pan has an amplitude of 
e~4x 10~ 7 and an initial wavenumber of \k \a ~ 1100, the surface density variations due 
to this wave are quite small, only Acr/00 — ea\k\ ~ 4 x 10~ 4 , which again is too small 
for detection. However, it might be easier to see this wave further downstream, due to the 
increase in \k\ with distance x. For instance, \k\ will have increased by a factor of 10 when 
the wave reaches the outer part of the A ring, so Aa/a ~4x 10 -3 there, but this again 
is probably still be too small for detection. Thus the observational consequences of these 
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density waves are seemingly slight. Nonetheless, these waves are not totally inconsequential, 
since Section [3] will show that the excitation of these waves can also result in a vigorous 
damping of the satellite's eccentricity. 



2.3.1. viscous damping of spiral density waves 

The variations in the ring's surface density due to the density wave are dominated by 
the second term in Eqn. (1531) . so Act ~ — ooefca sin(0 — uj), where the ring's longitude of 
periapse is u>(a,t) = out — f a k(r)dr, with the out term accounting for the spiral pattern's 
rotation with time due to the system's precession. Inserting this into the above shows that 
Act has the form Act = 3? (Se %< ^), which has an amplitude |£(a)| = cr e|/c|a and a phase 

0(a, 9, t) = mO - cut + J k(r)dr, (55) 

noting that 9 in the above was multiplied by m = 1 so that this work is readily compared 
to other studies of density waves developed for spiral patterns having m > 1 arms. 

When there is no dissipation in the system, the wavenumber k in Eqn. ( 1551) is real (e.g., 
Eqn. Wf\ . However, if there is any dissipation in the disk, then the wavenumber k acquires 
an imaginary component, i.e., k — > k^ + iki, which causes exponential damping of the wave's 
amplitude. If that dissipatio n is due to the disk's kinematic viscosity u, then the imaginary 
part of the wavenumber ia3 ( Shu 1984 ) 



vk\ 7skvk 2 R (mn — uj) vk R ( _ 7 



fc/ = -^+ K — — [ 1 + 7 — rr~ r J (56) 

mn — 00 D7tGcto n \ bjJd\kR\a J 

where kn is the real part of the wavenumber (Eqn. 14"?]) . = sgn(A^) = +1, m — 1, and 
noting that nGcr = fi d an 2 and that the spiral pattern rotates slowly, i.e., uj <^ n. Since 
the long waves launched by Pan have wavenumbers of \k\a ~ 10 3 to (see Section P2.3h . and 
that the A ring's normalized disk mass is fid ~ 4.5 x 10~ 8 (Section 12. 2. 4p . it is clear that the 
second term in Eqn. (|56|) dominates over the first, so 

kja ~ — (-?-) \k R a\ 2 (57) 



(5? 



6 fid \a?n 

for outward-propagating long waves. Consequently, the ring's surface density varies as 



Act = 3? (Se i(t> ) oc exp f- k^drj 



1 Comparison of Eqn. ([55)) to the phase convention adopted in IShul (jl984l ) shows that the signs of m and 
to are reversed. This sign reversal is accounted for in Eqn. (1561) . 
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which reveals how ring viscosity reduces the amplitude of the wave as it travels a fractional 
distance x. 

It is anticipated that the wave's viscous damping length l v will be much larger than 
the gap half-width A. In that case, the wavenumber that appears in the above formulae is 
approximately 



\kn\a 



(59) 



tt/^A 2 

by Eqn. (J37J). Inserting this into Eqn. (|57p and then evaluating the integral in Eqn. 
shows that 

Aa oc e - (x/c)3 (60) 
where i v is the viscous damping length in units of a s : 



a 2 n\ (ixDA 2 



18^| 
7 



The following evaluates this viscous damping length 
at the outer edge of the Encke gap. 



1/3 



(61) 



for the waves that Pan would launch 



There are two sources of viscosity in Saturn's rings: collisions among ring particles, and 
gravitational wakes. Viscosity due to collisions among ring particles is v c ~ 0.46c 2 r/n(l + 
t 2 ) where c = hn is the p article's dispersion velocity and r is the ring optical depth 



( jGoldreich fc Tremaindll982l ). Thus the dimensionless viscosity combination u c /a 2 n in the 
above evaluates to 

v c 0.46hV 



a?n 



x 10 
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(62) 



since the outer A ring has an optical depth of r ~ 0.6 and f) ~ 2 x 10 7 . Note that the 
viscosity that is associated with the ring's gravitational wakes v g is also comparable, since 



a?n 



CG 2 a 2 Cfi 2 



a 2 n 4 



7 x 10" 15 



7T- 



(63) 



where the coefficient C — 33 for the A ring ( Daisaka et al.ll200ll ) and /x^ ~ 4.5 x 10 -8 . This all 
suggests that the total viscosity in the outer A ring is around v/a 2 n ~ 10 -14 . Inserting this 
into Eqn. ( IBTT) then shows that the viscous damping length for Pan's waves (with A = 0.0012 
and /i c = 7.9 x 10~ 8 ) is t v ~ 0.039, which corresponds to a physical distance of £ u a s ~ 5200 
km. Note, though, that the outer edge of the A ring only lies a fractional distance distance 
of x ~ 0.024, which means that the amplitude of Pan's outbound density waves is reduced 



only by a factor of e 



0.8 as it travels across the A ring. This calculation indicates 



that the damping of long waves due to ring viscosity is only of marginal importance. 
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Section 14731 will also show that this long wave will reflect at the outer edge of the A ring 
(or else at the nearby Keeler gap) and march back towards Pan as a superposition of both 
long and short waves that are of roughly equal magnitudes. Pan's outbound long density 
waves will have a dimensionless wavenumber H L = y/2t)\k\a ~ 3.4 x 1CT 3 <C 1 when they 
near the outer edge of the A ring. Let H$ be the dimensionless wavenumber of the reflected 
short waves. Since the long and short spiral wave patterns both precess at the same rate, 
the wavenumbers Hl and Hs both satisfy the same dimensionless dispersion relation, Eqn. 
(I35D, so 

H L ~u* = H s e- Hs . (64) 

This is solved numerically for the wavenumber of the reflected short wave, Hs — 7.7, which 
yields a physical wavenumber of \k s \a ~ 2.7 x 10 7 since f) ~ 2 x 10~ 7 . This corresponds to a 
wavelength X s = 2ir/\k s \ that is comparable to the disk's scale height h = i)a s ~ 30m. These 
very short- wavelength waves will be nonlinear (Aer/er > 1), and they will likely damp on a 
very short spatial scale. 



3. Damping the satellite's eccentricity 

The excitation of these density waves also alters the satellite's eccentricity e s at a rate 
that can be calculated from the angular momentum flux that is transported by these waves. 
The angular momentum content of a narrow annulus in the disk is SL = 5m^GMa(l — e 2 ), 
where 5m = a5A is the mass of that annulus which has an area 5 A and a mass surface 
density a. The surface density of angular momentum in that annulus is then I = 5L/5A = 
ana 2 \/l — e 2 ~ ana 2 (l — e 2 /2) = £o + £ w , where Iq = ana 2 would be the angular momentum 
surface density if the ring were circular, and l w = —ae 2 na 2 /2 is the surface density of angular 
momentum that is associated with the wave whose amplitude is e. The flux of angular 
momentum that is transported by the wave is F = £ w c g , where c g = Hd<in is the group 
velocity of long density waves (see Section |2.2.2p . Consequently, the angular momentum 
luminosity, which is the rate at which waves transport angular momentum across an annulus 
of radius a, is £ = 2naF = — (e/^an) 2 M. 

Note the sign on C, which means that these waves transport angular momentum inwards, 
from the disk to the satellite. This transport also increases the satellite's angular momen- 
tum L s = m s yjGMa s {l — e 2 ) ~ m s n s a 2 (l — e 2 /2) at the rate dL s /dt ~ —m s n s a 2 e s e s = —C, 
which then provides the rate at which wave excitation tends to damp the satellite's eccen- 
tricity: 

p 2 n 2 

e s = --^n s . (65) 
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And if the wavenumber varies little over that first wavelength, which requires fi c <C fid, then 
the wave amplitude e is given by Eqn. fl52l) . which yields the satellite's eccentricity damping 
rate 

Vs(s + fic/ 'fid) 2 



( 1± 



and the eccentricity decay timescale 



(2ttDA) 2 
2tt(DA) 2 



where P or b 



fi s {s + fl c j 'fl d ) 2 
2n/n s is the satellite's orbit period. 



orb; 



(66) 



(67) 



If, however, fi c > fid, then Eqn. (152]) will overestimate the wave amplitude e, and Eqn. 
fl66l) will not be valid. Nonetheless, Eqn. (l65l) can still be used to determine the satellite's 
e-damping rate, but the wave amplitude e must be determined by other means, such as by 
using the rings model that is described in Section HI 



3.1. comparison to Lindblad and corotation resonances 

The satellite is also perturbing the planetary ring at its many Lindblad and corotation 
resonances in the disk, and this interaction also causes the satellite's eccentricity to evolve 
at the rate 

&s\lc = ^ s \l + (68) 
where the first term represents the eccentricity excitation that is due to the Lindblad reso- 
nances, and the second term is the eccentricity damping that is due to corotation resonances. 
The rate at which the satellite's eccentricity e s evolves due to its Lindblad resonances in a 
narrow annulus of mass 8m! that lies a fractional distance x' away is 

, _ f L ^e s W , , 

d6slL ~ 2|x'l 5 M [ } 



(from iGoldreich fc Tremaind Il98ll ). where the factor fi ~ 3.045. The satellite's total ec- 
centricity variation due to all of its Lindblad resonances throughout the disk is obtained by 
setting 5m' = 2-na'a'da' = 2fi' d Mdx' and integrating Eqn. ( |69l) across the disk, which yields 



CO 



es| L = 2 J a 5 ^\l = ^ L ^f S n s, (70) 

where the factor of two on the middle term accounts for the ring material orbiting interior 
and exterior to the satellite. The e-excitation timescale due to Lindblad resonances is then 

e A 4 

TL = -^r = —. Porb- (71) 
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The satellite also has many corotation resonances that lie in the planetary ring, and 
their effect is to damp the satellite's eccentricity at a rate e x \ c that has the same form 



as Eqns. flo^HTD]) but with a lead coefficient of fc = —3.193 (IGoldreich &: Tremaindll98ll ). 
However this eccentricity damping only occurs if the particles' motions at the corotation 
resonances are not saturated. The total rate at which the satellite's eccentricity varies due 
to its Lindblad and corotation resonances is 

e s\LC ~ e slL ' e s\C ~ 2 A 4 Us 

where the factor f LC = /l + fc — —0.148, which indicates that the net effect of these 
resonances is to damp the satellite's eccentricity, provided the corotation resonances are not 
saturated. Note, though, that \ fhcl fh\ — 5%, which means that eccentricity damping wins 
by only a small margin. Consequently, the total e-damping timescale due to Lindblad and 
corotation resonances is then 



tlc 



LC 



A 4 

^orb- (73) 



n\hc\lJ>sHd 



If, however, the corotation resonances are saturated, then the Lindbla d resonances will 



pump up the satellite's eccentricity, but at a rate that is half of Eqn. ( 1701) (see IGoldreich &: Tremaine 



19811 ). Nonetheless, the satellite's secular perturbations of the ring will also endeavor to damp 
the satellite's eccentricity at the rate given by Eqn. fl66l) . Comparing these two rates will 
show that eccentricity damping is still assured when the satellite's gap is sufficiently wide: 

A 2 > f L (7rD) 2 fi d /2, (74) 

provided // c <C A*d- 

Daphnis satisfies the /i c <C fid condition, and since fid — 4.5 x 10~ 8 , D = 0.87, and 
e = 2, the above gap-width requirement evaluates to A > 7 x 10~ 4 . However the Keeler 
gap's half width is only A ~ 1.5 x 10~ 4 , so this requirement is not satisfied. Consequently, 
if Daphnis' corotation resonances were in fact saturated, then this secular e-damping could 
not counterbalance the effects of the Lindblad resonances, which would pump up Daphnis' 
eccentricity until that satellite crashed into the A ring. That Daphnis has a very low eccen- 
tricity suggests that the motions of the ring particles at its corotation resonances are not 
saturated. 



Note that Pan does not satisfy fi c fid, so Eqn. f|74|) does not apply to this satellite. 
Instead, the rings model of Section 14731 is used to assess that satellite's e-damping rate. 
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4. Simulations of density waves launched at a gap edge 



The rings model of iHahnl (120031 ) is used to test the preceding results. This model treats 
the disk as a set of N discrete gravitating annuli having semimajor axes aj, eccentricities ej, 
longitudes of periapse u)j, and half-thicknesses hj. The model only considers the system's 
secular gravitational perturbations, so it also solves the same equations of motion, Eqns. 
(PQ), but it does so without making any of the approximations and assumptions invoked in 
Section [2j Consequently, the model provides another check on the analytic results obtained 
above. 



4.1. outbound waves 

The rings model is used to simulate the spiral density waves that are launched by an 
eccentric satellite that orbits just interior to a disk. Figure [3] shows the amplitude of this 
wave as it advances across a disk, with the system's parameters being detailed in the figure 
caption. Note, though, that those parameters do not correspond to any real ring-satellite 
system. Rather, these parameters were chosen to illustrate the results of Section [2] in the 
limit in which those results were obtained, namely, that the factor \x c j [id appearing in the 
wavenumber Eqn. (j!7|) is modest, with \i c j \i<i = 0.2. This causes the wavenumber to vary 
only a small amount across the first wavelength, which is a key assumption of Section 12. If s 
derivation of the wave amplitude. Nonetheless, the simulation reported in Fig. [3] does 
correspond loosely to a small ~ 10 km satellite orbiting just interior to a ring whose surface 
density is similar to Saturn's main A ring. 

The amplitude of the simulated wave compares favorably with the wave's expected 
amplitude, Eqn. (|52|) . which is indicated by the dashed line in Fig. [3j The wave propagation 
time also provides another check on these calculations. Note that the time for these waves 
to propagate a fractional radial distance x = Ar/a is 

c g 2vr/i d 

where c g is the waves' group velocity, Eqn. (|37|) . The simulated disk has a normalized mass 
of Hd = 5 x 10~ 8 and a fractional width x = 0.015, so the anticipated propagation time is 
tprop = 48 x 10 3 orbits, which compares well with with the simulated wave's propagation 
time (Fig. [3]). 

Section 12.2.21 predicts that the satellite will launch a long trailing density wave that has 
Sfc = sgn(fc) = +1. So by Eqn. f|T3|) . this means that the disk's longitude of periapse d)(a) 
should steadily decrease at greater distances a in the disk. This is confirmed in Fig. HI which 
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shows the waves' longitudes relative to the satellite's, Co {a) — Cb s . Also note that the longitude 
of periapse at the disk's inner edge is 90° behind the satellite's longitude, as expected. We 
find that once the density wave is established in the disk, its longitudes d)(a) precess at the 
same rate as the satellite's, and that the disk's eccentricities e(a) are also constant, which 
confirms the steady-state assumption employed in Eqns. The disk's surface density 
variations due to the wave is Eqn. (153]) . which is dominated by the second term. Inspection 
of that term shows that the disk's maximum surface density occurs at longitudes that trail 
those shown in Fig. H]by 90°, while the disk's minimum surface density Fig. H] leads by 90°. 

Figure H] also plots the dimensionless wavenumber |/c|oA across the disk at time t = 
60 x 10 3 , when the density wave just starts to reflect at the disk's outer edge. Also plotted is 
the expected wavenumber, Eqn. (|47|) . which compares well. The simulation's variations in k 
at the disk's far edge are due to the wave reflecting there, while the variations near the disk's 
inner edge are due to the very short wavelength variations in Co {a) that are just barely seen 
in the upper figure. Those wiggles are the short waves that are described in Section 12.2.21 
Evidently, the satellite also launches short-wavelength leading waves at the disk's inner edge, 
which are the cause for the high-frequency wiggles also seen in the t = 60 curve shown in 
Fig. [3J Figure H] also shows that the short waves propagate much more slowly than the long 
waves, as expected. Note, though, that nonlinear effects that are not modeled here would 
damp these nonlinear short waves on a very short spatial scale (e.g., Section [2.3. ip . 

The rate at which the disk damps the simulated satellite's eccentricity e s is shown in 
Fig. El which compares well with the expected rate, Eqn. 



The numerical quality of this simulation is assessed by monitoring the systems' total 
angular momentum deficit, which is L e = | ^\ rrijiija^e^, where the sum runs over all ring s 



and satellites in the system. Since this quantity is conserved by Eqns. (J2J) ( jHahnl 120031 ) . 
it provides a useful check on the simulation's numerical precision. The single-precision 
calculation shown in Fig. [3] conserves L e with a fractional error of \ AL e /L e \ < 4 x 10 -5 . 



4.2. reflection at an outer boundary 

Sect ion f2 .3.11 showed that ring viscosity is only marginally effective at damping the long 
waves that an eccentric satellite would launch at the gap's outer edge. Consequently, these 
waves will propagate outwards until they hit a barrier, such as another gap in the ring, or 
the ring's outer edge. The subsequent fate of such a wavetrain is illustrated in Fig. [61 which 
shows the state of the simulation described in Section H~T1 but at a later time, t — 2.1 x 10 5 
orbits. Here, the outbound long wave has already reflected at simulated ring's outer edge 
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and propagated back towards the satellite, but as the superposition of a long leading wave 
and a short trailing wave. Consequently, the long-wavelength undulations seen in Fig. [6fs 
e(a) curve represents the superposition of outbound and inbound long waves. Note, though, 
that the reflected waves' angular momentum content is now shared between a long and a 
short wave, so the amplitude of the inbound long wave is smaller than the outbound long 
wave. Consequently, the broad e(a) undulations seen in Fig. [6] are not due to a standing wave 
in the disk, but instead represent the superposition of two traveling waves having different 
amplitudes. 

The high-frequency variations seen in Fig. [6js e(a) curve are due to short waves that 
have a wavelength that is comparable to the disk scale height h. Note that the right-hand 
side of Fig. shows that the reflected short wave is still confined to the vicinity of the ring's 
outer edge by time t — 2.1 x 10 5 orbits, which indicates that these short waves travel slower 
than the long waves, as expected. The rapid variations in e(a) seen at the left side of Fig. [6] 
shows that the satellite is also exciting short waves at the gap's inner edge. These short waves 
are propagating much more slowly (see Section 12.2.21) . by a factor of He~ H ~ 0.01, so these 
very slow-moving waves have had little opportunity to travel very far by time t — 2.1 x 10 5 
orbits. Longer-term simulations also show that the amplitude of any short waves is always 
comparable to the amplitude of the long waves. Equation fl54"j) also indicates that these short 
waves will be very nonlinear (\Aa/a\ > 1), so the reflected short wave is expected to damp 
over a very short spatial scale that is probably comparable to the ring's scale height h, in 
the vicinity of the ring's outer edge. 

However, the reflected long wave will still return to the disk's inner edge, where it 
can interact with the satellite and/or reflect again. The satellite can interact with this 
returning wave by absorbing some of the wave's angular momentum content, which would 
also excite the satellite's eccentricity and seemingly stall any further e-damping. However 
that e-pumping would then be counterbalanced by enhanced eccentricity damping due to 
the excitation of even higher-amplitude density waves. Consequently, the damping of the 
satellite's eccentricity by this phenomenon is still assured, despite the fact that long waves 
propagate with only a modest amount of viscous damping. This is due to the fact that all 
long waves eventually reflect somewhere in the disk and spawn short waves that are easily 
damped. 

4.3. waves launched by Pan 

The rings code is also used to simulate the density waves that an eccentric Pan can 
launch at the outer edge of the Encke gap. As Section 12.31 notes, the ratio of the critical 
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mass fi c to the disk mass fid is nd H>d — 1-8, which indicates that the wavenumber changes 
substantially over the first wavelength (see Eqn. W7\ which violates a key assumption in the 
derivation of the wave amplitude. In fact, a comparison of Eqn. (|52|) to simulations of Pan's 
waves shows that that equation overestimates the wave amplitude by a factor of 7 ~ 4, which 
means that Eqn. (|66|) . which provides the satellite's eccentricity-damping rate, would also be 
in error by a factor of 7 2 . Nonetheless, a comparison of the simulated wave's wavenumber 
I A; I to Eqn. (T4T|) shows that equation to be in excellent agreement with the model results. 
Eqn. (I65p is also shown to be a reliable indicator of the satellite's eccentricity damping rate, 
even when /j lc > fi^. 

These same simulations also show that density waves launched by Pan are not able to 
propagate across the Keeler gap, which lies about 2900km downstream of the Encke gap. 
The width of that gap is approximately 40km, which is about half wavelength of Pan's waves 
in this region. Despite having a wavelength that is larger than the gap's full width, Pan's 
waves are unable to propagate across the Keeler gap. Instead, those waves reflect at the 
gap's inner edge, which, in the simulation, propagate back to Pan as a superposition of long 
and short waves. However in a real disk, those nonlinear short waves would quickly damp 
very near the Keeler gap's inner edge. 

This wave-action will also damp Pan's eccentricity over a timescale equal to Eqn. fIBTj) 
multiplied by r y 2 = 16, which evaluates to r e ~ 1400 years. Note that this e-damping 
is competitive with, but not quite faster than, the eccentricity excitation that is due to 
Pan's Lindblad resonances in the ring, which pump up the satellite's e over a slightly faster 
timescale of 77 ~ 900 years (see Eqn. [7T]) . But if the particles' motions at Pan's corotation 
resonances are saturated, then e-damping by th e corotation torque is shut o ff, while the e- 



excitation due to the Lindblad torque is halved (jGoldreich fc Tremaindll98ll ). so r L — > 1800 
years. So when Pan's corotation resonances are saturated, Pan's e-damping due to its secular 
interaction with the ring will exceed the e-excitation that is due to its Lindblad resonances in 
the ring, but only by a small margin. But if particle motions at Pan's corotation resonances 
are in fact unsaturated, then the near cancellation of the Lindblad and corotation torques 
results in a secular e-damping timescale that is about 13 times shorter than the total resonant 
e-damping timescale (Eqn. [731 . In this case, e-damping by the secular interaction is the 
dominant process that stabilizes Pan's eccentricity. 

Since this eccentricity damping due to wave excitation is so vigorous, one might w onder 



why Pan would even have a nonzero eccentricity. However, as ISpitale et al.l (120061 ) point 
out, the satellite Prometheus has a 16:15 resonance that does disturb Pan, which may be 
responsible for sustaining Pan's eccentricity and the density waves that that satellite would 
launch at the outer edge of the Encke gap. 
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5. Summary and Conclusions 

The Lagrange planetary equations are used to study the secular evolution of a small 
planetary satellite as it orbits within a narrow gap in a broad, self-gravitating planetary ring. 
These equations show that an eccentric satellite's secular perturbations of the nearby gap 
edge tend to excite very long-wavelength spiral density waves that propagate out to greater 
distances in the ring. These results are applied to the two small Saturnian satellite's Pan and 
Daphnis, which inhabit narrow gaps in the main A ring. It is shown that these satellites can 
launch very low-amplitude (Aa/a < 0.4%) long waves whose wavelengths would be of order 
A ~ 100's of km. The wavelength of these waves also shrinks with distance due to the central 
planet's oblateness, which causes the spiral pattern to wind up as the waves propagate. It is 
also shown that these long waves suffer only a modest amount of viscous damping as they 
propagate towards the A ring's outer edge. A dispersion relation is derived for these waves, 
which shows that a gap-embedded satellite can also excite short waves whose wavelength 
is comparable to the ring's scale height h. However these short waves are very nonlinear 
(Aa/a ~ 1), and will damp soon after their excitation. 

The Lagrange planetary equations are also used to derive the amplitude and wavelength 
of the long waves, as well as the rate at which this wave excitation tends to damp the 
satellite's eccentricity. However these analytic results are only valid when the wavenumber 
k varies slowly across the first wavelength, which requires the so-called critical mass \x c to be 
sufficiently small. Equation (j46p shows that that occurs when the gap is sufficiently narrow, 
or when the central planet's oblaten ess is sufficie ntly small. If, however, this requirement is 



not satisfied, then the ring's model of lHahnl (120031 ) can still be used to determine these waves 



properties, and the rate at which this wave action also damps the satellite's eccentricity. 
Note that the amplitude of these waves is proportional to the satellite's eccentricity, so this 
e-damping also tends to terminate subsequent wave generation. 

The rings model shows that these undamped long waves will eventually reflect at the 
ring's outer edge (or at another gap in the ring), which then spawns both long and short 
waves that propagates inwards. Since wave reflections eventually transmute all long waves 
into easily-damped short waves, this process communicates the wave's angular momentum 
content to the ring itself, which also insures that this wave phenomenon ultimately damps 
the satellite's eccentricity, too. 

This eccentricity damping due to wave excitation is then compared to the e-evolution 
rates that are due to the satellite's interaction with ring material orbiting at its Lindblad 
resonances (which tends to pump up the satellite's eccentricity) and corotation resonances 
(which tend to damp the satellite's eccentricity). It is shown that e-damping due to wave 
excitation is the dominant process when the gap width is sufficiently wide (see Eqn. I71"|) . For 
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the case of Pan, these e-damping and e-excitation rates are all comparable to each other, 
but Daphnis' e-evolution is dominated by the Lindblad and corotation torques. And since 
Daphnis' long-term orbital stability requires e-damping to dominate over e-excitation, these 
results also imply that particle motions at Daphnis' corotation resonances are unsaturated, 
which is necessary for the corotation torque to be operative here. 
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Fig. 1. — Schematic showing the geometry of the ring-satellite system seen edge-on. A 
satellite of mass m s and semimajor axis a s orbits interior to a broad planetary ring that 
extends to infinity. The satellite's distance from the ring's inner edge is A in units of the 
satellite's semimajor axis a s . A perturbed annulus in the ring has mass 5m, and it lies a 
fractional distance x away from the ring's inner edge, while the perturbing ring has mass 
5m! and lies a fractional distance x' from the perturbed ring. 
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I ■ ■ ■ ■ 

12 3 4 

wavenumber H 

Fig. 2. — The dimensionless dispersion relation, Eqn. (I35l) . plotted versus the dimensionless 
wavenumber H = y/2i)\k\a. Note that u>* has a maximum value of cu^ ax = exp(— 1) ~ 0.368, 
which occurs when the wavenumber H — 1. When the spiral density pattern precesses at an 
angular rate u* < u*^, this dispersion relation has two solution: a long-wavelength wave 
that has a wavenumber H < 1, and a short-wavelength wave that has wavenumber H > 1. 
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Fig. 3. — The rings model is used to simulate spiral density waves launched by an eccentric 
satellite that orbits just interior to a disk. The satellite's normalized mass is // s = 10~ 12 , 
and the disk is comprised of N = 500 rings having semimajor axes distributed over 1 + A < 
cij/a s < 1.015, where A = 5 x 10~ 4 is the fractional distance between the satellite and the 
innermost ring. The rings' fractional masses are \x r = 2.9 x 10~ 12 , so the normalized disk 
mass is fid = 7rar 2 /M = (fi r /2)(a s /5) = 5 x 10~ 8 , where the rings' fractional separations 
are S/a s = 0.0145/iV = 2.9 x 10~ 5 . The rings' fractional half-widths 1} is also set equal to 
their separations S/a s . The central planet's zonal harmonic is J2 = 0.012 and the planet's 
radius is R p /a s = 0.45, so this system's critical disk mass is \i c = 1.0 x 10~ 8 and \i c j = 0.2. 
The satellite's initial eccentricity is e s = 10~ 5 , with all other rings initially having zero 
eccentricities. The curves show the fractional amplitude of the density wave, e(a)/e s , as 
it advances across the disk, shown at selected times t in units of 10 3 orbital periods. The 
dashed line is the expected wave amplitude, Eqn. fl52l . with e — 1. 
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Fig. 4.— The upper figure shows the disk's longitude of periapse a) (a) relative to the 
satellite's periapse u s , in units of tt, for the simulation of Fig. [3] at time t = 60 x 10 3 orbits, 
when the wave has swept across the disk. The dots in the lower figure shows the dimensionless 
wavenumber |/c|aA at this moment, where wavenumber is calculated from k = —duo /da. The 
grey line is the expected wavenumber, Eqn. ( H7I) . with e = 1. 
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Fig. 5. — The rate at which the satellite launching the wave in Fig. [3] has its eccentricity 
damped, e s , plotted versus time t in units of 10 3 orbit periods. The solid gray curve is the 
expected rate, Eqn. ( |66l) . 
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Fig. 6. — The simulation described in Section 14.11 and Figs. [3HS] is extended to time t = 
2.1 x 10 5 orbits. During this time, the outbound long wave has since reflected at the simulated 
disk's outer edge at a/a s = 1.015 as a superposition of inbound long and short density 
waves. The long- wavelength undulations seen in e(a) thus represent the superposition of the 
outbound and inbound long waves, while the high-frequency variations in e(a) are due to 
the slower moving short waves. See Section 14.21 for details. 



